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Abstract: We study the spectrum of gravitational perturbations around a vacuum 
de Sitter brane in a 5D asymmetric braneworld model, with induced curvature on 
the brane. This generalises the stealth acceleration model proposed by Charmousis, 
Gregory and Padilla (CGP) which realises the Cardassian cosmology in which power 
law cosmic acceleration can be driven by ordinary matter. Whenever the bulk has 
infinite volume we find that there is always a perturbative ghost propagating on 
the de Sitter brane, in contrast to the Minkowski brane case analysed by CGP. We 
discuss the implication of this ghost for the stealth acceleration model, and identify 
a limiting case where the ghost decouples as the de Sitter curvature vanishes. 
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1. Introduction 

There is now a wealth of evidence [[IJ, [| indicating that the universe is currently 
in a period of accelerated expansion. One of the biggest challenges in cosmology 
today is understanding the origin of this late time acceleration. One possibility is 
that 70% of the energy content of the universe is dominated by an as yet unknown 
form of energy, so-called dark energy. The most popular dark energy candidate is the 
vacuum energy, which takes the form of a small and positive cosmo logical constant. 
In order to explain the current acceleration, the value of the cosmological constant 
must contribute a vacuum energy density of the order pa ~ lCT 12 (eK) 4 . This is 10 120 
times smaller than what we might expect, given our current understanding of particle 
physics. Given that particle physics is doing such a miserable job of explaining the 
accelerated expansion, it is important to look for alternative explanations. 

A popular alternative is to interpret this acceleration as a sign that our un- 
derstanding of gravity is breaking down, and that a large distance modification of 
Einstein's General Relativity is required. Despite numerous attempts, it is fair to 
say that an established proposal has yet to emerge that is consistent on both a fun- 
damental and a phenomenological level. Arguably the most successful attempts have 
been inspired by the braneworld paradigm (for a review see [[|). In particular, the 
Dvali-Gabadadze-Porrati (DGP) model || was discovered to have two cosmological 
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branches, one of which gave rise to cosmic acceleration even when no matter was 
present on the brane || . This branch became known as the self-accelerating branch, 
for obvious reasons, but was later discovered to be haunted by ghost instabilities 
around the vacuum de Sitter brane ||, [| |HJ H], [12|, [1^, |14], [15], |TB, [T7] (for a 



review see |n§). In this context a ghost is a field whose kinetic term has the "wrong" 
sign. This pathology leads to a choice: either the ghost state has negative norm 
and unitarity is violated, or the ghost can have arbitrarily negative energy. A ghost 
in the perturbation spectrum (spectroscopy!) indicates a catastrophic instability of 
the background, and therefore an unacceptably sick perturbative theory. In DGP, 
the other cosmological branch (the "normal" branch), is ghost-free but cannot be an 
alternative to ACDM since it still needs the introduction of the cosmological con- 
stant A to explain the acceleration. Nevertheless it still has plenty of interesting 



phenomenological features [19, EH, 21 



More recently, Charmousis, Gregory and Padilla (CGP)p^] presented a general- 
isation of the DGP model in which they allowed for bulk curvature and introduced 



some asymmetry across the brane |23|, |23|, |25], Ej], EJ, |28|] . This asymmetry could, in 
principle, apply to the bulk cosmological constant or even the bulk Planck scales, giv- 
ing rise to a rich variety of cosmologies. The authors focussed on those solutions that 
possessed asymptotically Minkowski branes, despite the presence of self-accelerating 
solutions that they (correctly) assumed to be haunted by ghosts. A subset of these 
solutions were shown to contain vacuum branes that were perturbatively stable, free 
from the ghoulish instabitities that terrorized the self accelerating DGP brane. The 
cosmological evolution of this subset was then analysed, and in some cases yielded 
extremely interesting results. Two limiting models in particular (the "decoupled" 
limit and the "conformal" limit) were found to exhibit power law acceleration but 
only when matter is present on the brane. They dubbed this 'stealth acceleration'. 
The cosmology is reminiscent of the Cardassian cosmology proposed by Freese 



and Lewis ||29|| . Here the standard Friedmann equation is modified so that p — > 
p + cp n , where n < 2/3, and one also finds that cosmic acceleration is driven by 
the presence of ordinary matter. The Cardassian model is an interesting empirical 
model, but did not have a concrete theoretical basis. The stealth model provides that 
by realising an effective Cardassian cosmology (with n « 0.5) within the braneworld 
paradigm. 

In this paper we will consider vacuum de Sitter branes within the CGP set-up. 
This will include self-accelerating solutions, as well as the stealth models with some 
additional vacuum energy on the brane. We will study the spectrum of linearised 
perturbations about these solutions, closely following the corresponding analysis in 
the DGP model P, 1C, 13fl. For an infinite volume bulk, we will find, without ex- 
ception, that the vacuum is unstable because of the presence of ghosts. Just as for 
the self-accelerating branch of DGP, a ghost will manifest itself either through the 
radion mode, or through the helicity mode of the lightest graviton. In some cases 
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a ghost will also appear in the spin 1 sector. 

The "decoupled" version of the stealth model is now of particular interest. We 
will find a class of de Sitter solutions that approach the "decoupled" model as the 
Hubble scale H — > 0. As the limit is approached the ghost becomes more and 
more weakly coupled, until eventually it decouples completely. We will infer some 
conclusions regarding the stability of the stealth models when matter is present. For 
small H it seems that we can carry our analysis of de Sitter branes over to the 
general Friedmann-Robertson- Walker case and conclude that the decoupled stealth 
model develops an instability albeit a very mild one softened by the weakness of 
the ghost coupling. For larger H the instability for de Sitter branes would be more 
severe, but it is not clear whether or not we can transfer this conclusion to the general 
FRW case. 

The rest of this paper is organised as follows: in section Q we describe the 
CGP model in detail, our generalisation, and the background solutions. In section 
|3] we analyse the spectrum of linearised perturbations and derive conditions for the 
presence of an helicity ghost in the spin 2 sector. We study the coupling to matter in 
section |] and calculate the effective action in section |5|. The effective action helps to 
reveal any further ghosts, including the radion ghost, which seems to take it in turns 
with the helicity mode to haunt the background. We end with some concluding 
remarks in section ^| 

2. The CGP model: set up and background solutions 

The CGP model is an asymmetric generalisation of its celebrated cousin, the DGP 
model. In both models, our Universe is taken to be a 3-brane, S, embedded in 
between two five dimensional spacetimes, A^, where i = L,R. In the original DGP 
scenario, we impose Z 2 symmetry across the brane, identifying M. L with M.r and 
having vanishing vacuum energy in the bulk. In the CGP model, however, we relax 
both of these assumptions. The key new ingredient is the introduction of asymmetry. 
Each spacetime Aii generically has a five dimensional Planck scale given by Mj, and 
a negative (or zero) cosmological constant given by Aj = —Qkf,. However, since 
we are no longer assuming Z 2 symmetry across the brane, we can have Ml 7^ M# 
and Al 7^ A^j. Allowing for 7^ Ar is familar enough in domain wall scenarios 
J5T]]. The Planck scale asymmetry is less familiar, but could arise in a number of 
ways. Suppose, for example, that this scenario is derived from a fundamental higher 
dimensional theory. This theory could contain a dilaton field that is stabilised in 
different fundamental vacua on either side of S. From the point of view of a 5D 
effective description, the 5D Planck scales would then differ accordingly. Indeed 
naive expectations from string theory point towards this asymmetric scenario as 
opposed to a symmetric one. Different effective Planck scales can also appear on 
either side of a domain wall that is bound to a five- dimensional braneworld |J3|] . 
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In keeping with the braneworld paradigm, all matter and standard model in- 
teractions are confined to the brane, although gravity can propagate into the fifth 
dimension. As in the DGP scenario, we include some intrinsic curvature induced 
on the brane. This term is rather natural and can be induced by matter loop cor- 



rections p3[, finite width effects |34| or even classically from higher dimensional 



modifications of General relativity We will also include some vacuum energy 
on the brane in the form of some brane tension, a . At this point we introduce 
an important new development. In the original CGP paper, the brane tension was 
fine-tuned against the bulk cosmological constants in order to admit a Minkowski 
vacuum solution. This choice corresponds to having vanishing effective cosmological 
constant on the brane and was the analogue of the Randall-Sundrum fine-tuning. In 
this paper we will introduce some additional tension so that the vacuum brane is 
de Sitter. Such detuning of brane tensions helped conjure up the ghost in the DGP 
model, and we will ultimately find that the same is true here. 

This set-up is described by the following action, 
S= V Mf [ ^(R-2A l ) + 2M? [ y/=jK® + [ ^(MtK-a+jC m&ttei ), 

i=LR J Mi JdM t JT, 

(2.1) 

where g a b is the bulk metric with corresponding Ricci tensor, R. The metric induced 
on the brane is given by = g a t — n a nb where n a is the unit normal to dM.i in Aii 
pointing out of Aii. Of course, continuity of the metric at the brane requires that 
'jab is the same, whether it is calculated from the left, or from the right of the brane. 
In contrast, the extrinsic curvature of the brane can jump from right to left. In Aii, 
it is defined as 

K-ab — jllffl(c n d), (2.2) 



with its trace appearing in the Gibbons-Hawking boundary term in ( [2.1| ). In the 
brane part of the action we have included the brane tension, a, and the induced 
intrinsic curvature term, 71, weighted by a AD mass scale, M 4 . £ mat ter includes any 
additional matter excitations. 

The equations of motion in the bulk region, are just the Einstein equations, 
with the appropriate cosmological constant, Aj. 

E a b = Rab — -jRQab + Aj^afe = 0. (2.3) 

The equations of motion on the brane are described by the Israel junction conditions, 
and can be obtained by varying the action ( |2.1| ), with respect to the brane metric, 
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7afe- This gives 1 

e ab = 2 (M 3 (K ab - K lab )) + Ml hiab - \n la A + a - lab = ^T ab , (2.4) 

where T ab = — -/= ^~2^™ ttcr ■ Note that the Israel equations here do not use the 
familiar "difference", because we have defined the unit normal as pointing out of 
Aii on each side. We adopt this (slightly) unconventional approach since it is more 
convenient in the asymmetric scenario where the brane is best thought of as the 
common boundary £ = 8M.L = 9M.R. 

We will now derive the vacuum solutions to the equations of motion (|2.3|) and 
( p.4|) . This corresponds to the case where there are no matter excitations, and so, 
T ab = 0. In each region of the bulk, we introduce coordinates x a = (x M ,y), with the 
brane located at y — 0. We are interested in de Sitter brane solutions of the form 

ds 2 = g ab dx a dx b = dy 2 + N(y) 2 ^^dx^dx v '. (2.5) 

where 7 M „ is the four dimensional de Sitter metric with curvature, H. Inserting this 
into the bulk equations of motion (Ot) gives 



^V-^ + P — — k 2 (26) 

where "prime" denotes differentiation with respect to y, and we have dropped the 
index i for brevity. One can easily show that 

H 1 
N(y) = — sinh k (y h + 6y), y h = - sinh" 1 k/H, (2.7) 

where 9 — ±1. Each region of the bulk corresponds to < y < y max where 

{oo for 6 = 1, . 
* a 1 (2 ' 8) 
y h for 6 = -1. 

If we transformed to global coordinates in the bulk, 6 = 1 would correspond to 
retaining the asymptotic region (large radius), whereas 9 = — 1 would correspond to 
retaining the central region (small radius). For k ^ 0, this means that when 8 = 1 
we keep the adS boundary (growing warp factor) whereas when 9 = — 1 we keep 
the adS horizon (decaying warp factor). Since we are interested in a modification of 
gravitational physics in the infra-red, we will assume that the bulk volume is infinite, 
and retain the asymptotic region on at least one side of the bulk. In other words, we 
do not consider the case 9l = Or = —1- 



1 Thc angled brackets denote an averaged quantity at the brane. More precisely, for some quantity 
Qi defined on the brane in dAii, we define the average (Q) = Ql + Qr . Later on we will also make 
use of the difference, AQ = Ql — Qr. 
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The boundary conditions at the brane ( |2.4|) yield 



6(M 3 N'(0)) + | - m 2 M\ = 0, (2.9) 



so that the curvature H is given by the real roots of 



a = QMlH 2 - 12 ( M 3 6VH 2 + k 2 ) . (2.10) 



In 1 22fl ) the brane tension was fine tuned to a critical value, a c = —6(M 3 k), so that 
the effective cosmological constant on the brane vanished. We now introduce some 
additional tension e > so that a = o c + e. This introduces some positive curvature 



given by the roots of e = F(H 2 ) where, as in [22]|, we have 

F(H 2 ) = 6M%H 2 - 12 (m 3 6 U/H 2 + k 2 - k) } . (2.11) 

As in DGP, we have two classes of solution. There are those that vanish as e — > 0, 
so that we recover the Minkowksi brane studied in Iffifl , and there are those that 
approach a finite positive value, so that we have a de Sitter brane, even in the 
absence of an effective cosmological constant. The former are the analogue of the 
normal branch in DGP, whereas the latter are the analogue of the self-accelerating 
branch. Of course, the class of solution depends on the form of the function F(H 2 ), 
discussed in some detail in section 4 of fl22| . For example, the following represent 
necessary and sufficient conditions for the existence of a normal branch solution: 

Ml > (M 3 8/k), (2.12) 
or M 2 A = (M 3 6/k), (M 3 6/k 3 ) > 0, (2.13) 
or Ml = (M 3 9/k), (M 3 9/k 3 ) = 0, (M 3 6/k 5 ) < 0. (2.14) 

Although we will study both classes of solution, we will be particularly interested in 
the normal branch since these will include small fluctuations about the finely tuned 



"stealth" scenarios discussed in [22 



3. Vacuum fluctuations 

We shall now consider metric perturbations in the vacuum so that g a b = g a b + Sg a b 
and T^y = 0. In the unperturbed spacetime, given by ( |2.5| ) and (|2.9| ), the gauge 
was fixed in both M.\ and M.2 so that the brane was at y — 0. However, a general 
perturbation of the system must also allow the brane position to flutter. In Aii, the 
brane will be located at 

y = &(*"). (3.1) 
It is convenient to work in a Gaussian Normal (GN) gauge, so that in M. i we have 

s 9yy = Sg m = 0, 5g lu , = h illl/ (x,y). (3.2) 
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In most of this discussion, we will drop the index i although its should be understood 
that it is really there. Now, it is well known (see, for example, ||36|| ) that in the absence 
of any bulk matter, we may take to be transverse-tracefree D^h^ — M = 0. This 
is known as Randall-Sundrum gauge. It follows that the bulk equations of motion, 
SE ab = give 

d 2 + _L( D 2 



N 2 



AH' 



Ak 1 



Ku(x,y) = o, 



(3.3) 



where is the covariant derivative on the AD de Sitter slicings, and indices are 
raised/lowered using the AD metric j^. To impose the boundary conditions at the 
brane, we need to apply a GN to GN gauge transformation that shifts the brane 
position back to y — 0. The most general such transformation is given by 

" y dz 



x h 



x r 



Jo 



N 2 (z) 



so that 



h 



+ hf u +2N 2 D { ^ u) . 



h, 



(3.4) 



(3.5) 



We call this new gauge "brane-GN" gauge. Although the brane position is fixed in 
this gauge, the original position ((x) still enters the dynamics through a bookkeeping 
term 

h$ = -2 (V jT ^\ D,D U C + 2NN% V (. (3.6) 

The metric perturbation in the new gauge is no longer transverse-tracefree, although 
it is now straightforward to apply continuity of the metric at the brane 



AV(i,0) =0, 



and the vacuum Israel equations 



(3.7) 



50 



liv 



- ( M z 



u liv 



in' 



N 2 



y=o 



+ MlX^h) = 0, 



(3. 



where 



3H 2 h 



[IV 



■-{D 2 - 2H 2 )h^ + D^D a h u)a - -D^DJi 



If we substitute the expression 

'h 



--% v [D a D l3 h a f 3 -(D 2 + H 2 )h\ 
into equation ( |3.8|) we find 



iw 



ML 9 

y=0 2 



2H 2 )h flu (x,0) 



2(D,D U - (D 2 + 3tf 2 )V) <(M 3 - M 2 N'(0))() . 



(3.9) 



(3.10) 
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Note that this expression is independent of ^(x), as expected, since this just corre- 
sponds to diffeomorphism invariance along the brane. It is convenient to decompose 
h^ v in terms of the irreducible representations of the AD de Sitter diffeomorphism 
group 

hr = h$+h$+h<®, (3.11) 

where h$ corresponds to the spin n contribution. We can treat these modes inde- 
pendently of one another provided they have different masses 2 . Let us now assume 
that this is indeed the case and analyse each spin separately. It will also be convenient 
to decompose the field ^(x) into its spin 1 and spin components £ M = ffi + . 
The field ((x) is just spin 0. 

3.1 Spin 2 modes 

We begin by analysing the spin 2 modes. Since neither £ nor £ M have a spin 2 
contribution, we can set them zero here, and can further decompose the spin 2 piece 
of the metric by separating variables 

h$(x,y) = [ u m (y)xt\x), (3-12) 



where xffi is a AD tensor field of mass m satisfying (D 2 — 2H 2 ) x ffi (x) = m 2 x^(x), 
and J" denotes a generalised sum, summing over discrete modes and integrating over 
continuum modes. The bulk equations of motion (|3.3j) now give 



'Ou)+[ m \^ H2 -4^)»,JU)-- 0. (3.13) 



AT2 

This is easily solved in terms of the associated Legendre functions: 

/ i \ 2 / , \ 2 

u m (y) = Ci Ijj) Vtl /2±v (coth k(y h + By)) + C 2 l-\ Q ±2 l/2±u (coth k(y h + By)) , 

(3.14) 

where v = ^/9/A-m 2 /H 2 . V™{z) and (z) are the associated Legendre functions 
of the first and second kind, respectively. Of course, the expression ( |3.14j ) is only 



well defined for m 2 < We could, in principle analytically continue our solution 
to m 2 > although this will not be necessary since our ultimate goal is to 

establish the existence of an helicity-0 ghost which is found in spin 2 modes of mass 
< m 2 < 2H 2 |g . Normalisability requires that f2"7j 



£/max y2 







dy ^<oo, (3.15) 



so that for B = 1 we only keep the part proportional to T > _i/2+ v i. z )i whereas for 
B = — 1 we only keep the part proportional to Q 2 _ 1 / 2 + l/ ( z )- Since we may assume 



2 In AD de Sitter, a transverse-tracefree tensor of mass m satisfies (D 2 — 2H 2 )qffi = m 2 qffi 



-8- 



that u m (0) 
given by 



1, without loss of generality, we get that the normalizable modes are 

,(cothfc (y h +y)) 



u m {y) 



-1/2+1/ 



(coth k (y h -y)) 



for 9 
for 9 



+1, 
-1. 



(3.16) 



It will be instructive to take a closer look at two special cases. For massless modes, 
this expression simplifies to give 



My) 




2+coth k{y h +y) 
2+coth k yh 



N 2 /^ max dz/N 4 
/ ymax dz/N 4 



for 9 
for 9 



+1, 
-1. 



(3.17) 



whereas for "partially massless" modes of mass m? 

-2ky for q 



U 



V2H 



(y) 



NN' 
N'(0) 



for 9 



2H 2 we have 

+1, 
-1. 



(3.18) 



Of course, neither the massless modes, nor the partially massless modes get excited 
in general. This is determined by the boundary conditions at the brane. The spin 
2 part of the continuity equation ( |3.7| ) now implies that Axffife) = for each m, 
so that the spin 2 part of Israel equations ( |3.10| ) yield the following quantization 
condition 

^_ / M 3 



f(m 2 



\N 2 J 



y=o 



Ml 2 

H m 

2 



0. 



(3.19) 



Let us consider the lightest mode. For a finite volume bulk (0l — 9r — —1), it is 
well known that this mode is massless so that gravity looks four dimensional out to 
arbitrarily large distances. We do not consider this case here, and assume, without 
further loss of generality, that Or = +1. The lightest mode is now guaranteed to 
be massive. If the mass lies in the forbidden region < m? < 2H 2 , then this 
mode contains an helicity-0 ghost [30|. We can now check if such a mode exists, by 
application of Bolzano's theorem: 



f(0)f(2H 2 ) < 0, 



(3.20) 



since f(m 2 ) is continuous over the forbidden region. Although not necessary for 
the existence of a ghost, this condition is certainly sufficient. For an infinite bulk 
9r) 7^ (—1, —1)), it is easy enough to see that 



/(0) = - -(M 3 (l + 



Smax dz 



< o. 



(3.21) 



This means we have an helicity-0 ghost whenever 



f{2H 2 ) = — 



M 3 
2 



(1-9)H 2 
VH 2 + k 2 



2(1 + 9)(k + y/H 2 + k 2 ] 



+ MlH 2 > 0. (3.22) 
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3.2 Spin 1 modes 



We now turn our attention to the spin 1 modes, neglecting all contributions from 
spin 2 and spin 0. Recall that £ M contains a spin 1 piece £,jp(x), which is simply 
a divergence-free vector that can be chosen in order to guarantee continuity at the 
brane. The spin 1 part of the metric takes the form 



D fl A u 



(3.23) 



where A^ix^y) is another divergence free vector. Since h^} is transverse-tracefree, 
one can easily verify that A^ behaves like a tachyonic vector in dS^, satisfying 



(D 2 + 3H 2 )A^ 



0. 



(3.24) 



This tachyonic instability is a mild one, associated with the repulsive nature of 
inflating domain walls f3"5R . The metric contribution now resembles a massless spin 
2 mode, (D 2 — 2H 2 )h^l = 0, and is therefore guaranteed not to mix with any of the 
genuine spin 2 modes discussed in the previous section. Furthermore, it follows that 
the profile in the bulk is given by the normalisable massless wavefunction ( |3 . 1 7| ) 



Ap(x,y) = u {y)a p (x). 



(3.25) 



The spin 1 part of the continuity equation (|3.7| ), A(a M + £^) = , is trivially satisfied 
by choosing (x) 



are independent of ^ , and require that 



-a^(x) on both sides of the brane. The Israel equations ( |3.10| ) 



y=0 



a M (x) 



0. 



(3.26) 



If we assume, without any great justification, that Aa M = 0, if follows from ( |3.26| ) 
that f(0)a^(x) = 0, and so a M (x) = 0. However, in a generalised asymmetric scenario 
there is no reason to assume that the spin 1 mode is symmetric. More generally we 
can show that 



/(0)<a M > 



-A M 3 (l + 



'//max 



LJO 



dz 



-V 



Aa 



(3.27) 



which indicates that one spin 1 degree of freedom can, in principle, remain. 



3.3 Spin modes 

We conclude this section with a study of the spin modes, neglecting all contributions 
from higher spin. The brane bending piece ( now plays a role, along with the spin 
component of which takes the form ^ = D^i/j, where ip(x) will be chosen in order 
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to guarantee continuity at the brane. The spin part of the metric perturbation can 
be written in terms of a pair of scalars, y) and h^°'(x,y), like so 



4^ 



(3.28) 



From the transverse-tracefree property of hffi, it follows immediately that = 
and 

(D 2 + 4H 2 )$ = 0. (3.29) 

Again, we have a mild tachyonic instability associated with inflating domain walls. 
The metric contribution now resembles a "partially massless" spin 2 mode, (D 2 — 



2H 2 )h$ = 2H 2 h^J, which could, in principle, mix with one of the genuine spin 2 
modes discussed in section |3.1| . We will discuss this in more detail later on. Assuming 
for the moment that there is no issue with mixing, we conclude that the scalar's profile 
in the bulk is given by the partially massless wavefunction 



${x,y) = u^ H (y)(j){x). 



(3.30) 



The spin part of the continuity equation is split into a pure gauge part, and a 
conformally de Sitter part. Requiring continuity of both parts separately implies 
that 

A(0 + 2^) = O, A(# 2 + 2jV'(O)C) = 0. (3.31) 

The first condition can be trivially satisfied if we chose ip(x) = —<f>(x)/2. The Israel 
equations (|3.10p are independent of ip, and require that 



M 3 



U V2H \' 

N 2 ) 



y=0 



+ MjH 2 <j){x) ) = 2((M 3 - M 2 N'(0))(). 



(3.32) 



It follows from ( fOID and (|3T32|) that 



A0 = -—A \ecvip + k? 



a 



9(VH 2 + k 2 } + PA \ 6(VH 2 + fc 2 J , 

(3.33) 



where 







f(2H 2 



M 3 9 



(vWTk 2 



Ma 



AH 2 f(2H 2 



■A 



(k + VH 2 + k 2 f 
VH 2 + k 2 



(3.34) 
(3.35) 



Here we see that the fluctuation in the brane position sources the bulk mode 4>(x). 
We therefore associate it with the radion. Again, there is no reason to assume 
A(j) = 0, so that in general the boundary conditions leave us with up to two spin 
degrees of freedom. Note that both a and j3 diverge as f(2H 2 ) — > 0. This singular 
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limit corresponds to the case where there exists a genuine spin 2 mode with mass 
m 2 = 2H 2 . The divergence in a and (3 reflects the fact that the lightest spin 2 mode is 
no longer orthogonal to the spin contribution, and cannot be treated independently. 
The two modes mix and a more careful analysis is required. Finally, we also note 
that we can write a = —F'( y H 2 )/3f(2H 2 ), where F(H 2 ) is given by equation ( |2.11| ). 



4. Coupling to matter 



When we introduce some additional energy-momentum, T^ u , on the brane, the homo- 
geneous solution discussed in the previous section picks up an additional contribution 
that describes the responses of fields to the source on the brane, 



hu„{x,y) -> h^ v (x,y) +7r lxu (x,y), Q(x) -> Q(x) + ix l {x), 



(4.1) 



where ir^ u is transverse-tracefree. In analogy with the theory of ordinary differential 
equations, it is useful to think of the homogeneous pieces, h^ u (x,y) and d(x), as 
the "complementary functions" and the inhomogeneous pieces, 7i^ u (x,y) and n^x), 
as the "particular integrals". The "particular integrals" must be solutions to the 
following 



d 2 y + ^- 2 (D 2 - AH 2 ) - Ak 2 



Ku„(x,y) = 0, 



A [tt^(x, 0) + 2N'(0)j tlv n(x)] = 0, 



(4.2) 
(4.3) 



M 2 

+ J 1±(D 2 - 2H 2 ) tt^(x, 0) 
y=o 2 



2{D tl D v - (D 2 + 3iJ 2 )V) ((M 3 - M 2 N'(0))tc(x)) + -T, u . 

Tracing the two boundary conditions, and carrying out a little algebra, gives 

T 



A [N'(0)n(x)] = 0, (D + AH 2 ) (N'(0)n(x)) 



(4.4) 



(4.5) 



2F'(H 2 Y 

where F(H 2 ) is given by equation fl2.11|) . This completely specifies the 7Tj(x), since 
any homogeneous brane bending is already accounted for in the 0(x). We now turn 
our attention to the ir^. The traceless part of (|4.3j ) demonstrates that Air^x, 0) = 
0, whereas the Israel equation ( f4.4|) may be rewritten like so 

1 , , 



y=0 I 



(4.6) 



where is a gauge invariant brane stress energy perturbation defined as [|T^] 



3 

+ ^{D^D U - (D 2 + 3H 2 )^^) 



T 



D 2 + AH 2 



(4.7) 
(4.8) 
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It turns out that in Ml* 



7T 



Ml 



(x,y) = / dVV=f G^{x,y;x',O)r a0 {x'), 



(4.9) 



where G^l a ^(x } y; x', 0) is the relevant Green's function, satisfying 



St + - 4# 2 ) - ^ 



G^(x, y] x',0) = 0, 
A[G^(x,0;x',0)} =0. 



G^(x,y;x',0)\ 



N 2 



M 2 

+ -±(D 2 -2H 2 )G^(x,0;x\0) 

y=0 2 



(4.10) 
(4.11) 
8^{x-x') 
"(4.12) 



The Green's function can be expressed in terms of the wavefunctions u m (y) discussed 
in section ( |3.1|) . Defining the normalised wavefunctions Um(y) = M m Um(y), where 
the M m are chosen so that 



Mlu m {0)u n {V) + 2M 3 / dy ^ 



o 



-nn for discrete modes, 

5(m — n) for continuum modes, 

(4.13) 



we have 



G^(x,y;x',0) = - / x « (aOx*«°V) 



^(y)^(0) 

^2 _ m 2 



(4.14) 



Note that (L> 2 - 2H 2 ) X ( $ = p 2 X$, and x* Q/3 satisfies 



For more details on this construction, at least for DGP gravity, see section 3.3 of [13 . 



5. The effective action 

We now compute the effective 4D action of normalisable vacuum perturbations. This 
will enable us to identify any ghosts: pathological modes with negative kinetic terms. 
Of course, we already know that whenever f(2H 2 ) > a ghost haunts the helicity- 
sector of the lightest spin 2 mode. Our effective action calculation will reveal a 
generic spin-0 "radion" ghost in the opposite regime, ie when f(2H 2 ) < 0. 

We begin our calculation in bulk Randall- Sundrum gauge, so that the brane is 
positioned at y = ((x) and the metric perturbation is given by 

M*>J/)= / ^Uy)x^\x)+u (y)h^(x)+u^ H (y)h^(x), (5.1) 
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where 

b${x) = D^a u + D v a^ h$(x) = (D^D U + H 2 % v )(j>. (5.2) 

In computing the action, it is important to leave the 4D fields off-shell. In other 
words, we do not assume (D 2 — 2H 2 )xffl = m 2 Xi%\ (D 2 + 3# 2 )a M = 0, or (D 2 + 
4H 2 )<f) = 0. These equations should follow from variation of the action at the end of 
the calculation. 

Randall-Sundrum gauge is the correct gauge choice far from the brane, since it 
contains no pure gauge modes with a non-normalisable profile in the bulk. However, 
in order to compute the effective action, it is convenient to be in brane-GN gauge 
close to the brane so that it lies at y = and the 4D coordinates match on either 
side. This can be achieved whilst maintaining Randall-Sundrum gauge far from the 
brane, but only at a price: we are no longer everywhere Gaussian-Normal. We can 
transform to this "fixed wall" gauge from everywhere Randall-Sundrum gauge by the 
following gauge transformation 

y^y-v y (x,y), x 11 -> x 11 -^(x^y), (5.3) 

where 

V v (x,y) = < r)>*{x,y) = \ {z > 

[0 for y>y*, [0 for y > 

(5.4) 

where < < y max is some appropriately chosen finite distance. It follows that 

$9ab $9ab + 2V( a ^ 6 ), (5.5) 

where V is the covariant derivative for g^. This new gauge interpolates between 
Randall-Sundrum gauge deep inside the bulk and brane-GN gauge near the brane. 
As result, the metric perturbation along the brane is the same as in brane-GN gauge, 
with 

<h> = M*> °) + 2jV '(°) VC + 2D^ V) . (5.6) 
We now perturb the action to quadratic order 

6S=(m 3 J d 5 x^6g ab 6E ab ^ + i J d^x^S^SQ^, (5.7) 



where 5E ab and 50^ are the linearised bulk equation of motion (2.3) and vac- 
uum Israel equation Q2.4|) , respectively. Using ( |5.5|) , ( |5.6[) and the Bianchi identity 
W a 5E ab = 0, we find that 

5S = I d 4 x^/^/5C, (5.8) 
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where 



5C 



-M i 



|/max 



dy h^{x,y)5E flv {h) 



+ 2M i V a (x,0)6E ay (h) 



y=0 



- (h" u (x, 0) + 2N'(0)r u ( + 2D^ U) ) SQ^. 



(5.9) 



We cannot assume D^hff = D^hfjfJ = = 0, since these imply the on-shell 
equations of motion for a M and <fi. We therefore need the following expressions for 
SE ab and SQ^ U for a generic GN perturbation. 



5E, y (h) 

5Eyy(h) 

5Q 



N 2 ' 



:<9„ 



H 2 

W 2 



2k' 



N 2 



2N y 



N 2 



2iV 4 



(D^D V - (D 2 + 3H 2 )^ v )h 



(5.10) 
(5.11) 
(5.12) 



M 3 a, 



N 2 



- MiX^h) 



y=0 



+2(D^D U - (D 2 + 3H 2 )^ U ) <(M 3 - M 2 iV'(0)K) . (5.13) 

Making use of equations (ETI3D , (jOSD , (gg§ , (ggg) , flpg) , as well as the orthogo- 
nality condition 

, U m {y)u n {y) 



2M d 



dy- 



N 2 (y) 



+ Miu m (0)u n (0) 



0. 



m ^ n, 



we arrive at the following 4D effective Lagrangian 

5C = 5C 2 + 8£i + 8C , 

where the spin 2, spin 1 and spin contributions are respectively given by 
1 



(5.14) 



(5.15) 



SC 2 
SC 



dy^l + M 2 



o 



N\y) 



2 

1/ I 

4\M3(u /N 2 )>\ y=0 

9f(2H 2 



X {m) ^(D 2 -2H 2 -m 2 )x^\ (5.16) 



-i 



Aa^(D 2 + 3H 2 )Aa 



ft! 



(5.17) 



F'(# 2 



-(7) 



(0) 



A7A0 
4(7) 



(L> 2 + AH 2 



<0> + 



A7A0 
4(7) 



3# 2 /l 



-1 



and 



7 



7 = M ; 



A0( J D 2 + 4# 2 ) A0, 

1 + (k + ^/WT¥) 2 
Vtf 2 + /c 2 



(wo/iV 2 )'| y= o = - 



2 

1 + 



> 0, 



dz 



< 0. 



(5.18) 

(5.19) 
(5.20) 
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As we stated earlier, we do not consider a finite volume bulk (9l = Or = — 1). Let 
us now analyse the alternatives. We see immediately that there is a spin 1 ghost 
whenever Ol = Or = +1. There are two spin modes, roughly corresponding to 
the average radion, and the difference. The latter is never a ghost, whereas the 
kinetic term for the average radion is determined by the sign of f(2H 2 )/F'(H 2 ). 
Recall that for a well behaved cosmology we require that F'(H 2 ) > [J22| . For finite 
F'(H 2 ) > 0, it follows that we have a radion ghost whenever f(2H 2 ) < 0. In section 
|3.1| , we found that the lightest spin 2 mode contains an helicity-0 ghost in precisely 
the opposite regime, ie when f(2H 2 ) > 0. This is exactly the sort of behaviour found 
on the self-accelerating branch of DGP: a well behaved spin 2 sector corresponds to 



a pathological radion, and vice versa || [H| [13], [14 1 . 

When 0l = —1, Or = +1, the spin 1 mode, and the radion difference decouple 
completely. In contrast, the average radion typically remains in the spectrum, and 
we can draw similar conclusions regarding its stability as discussed in the previous 
paragraph for 0l = Or = +1. However, there are a few exceptional cases. F'(H 2 ) = 
and F'(H 2 ) — > oo ultimately correspond to the "stealth" scenarios identified in [fZ2|1. 
where the brane is Minkowski as opposed to de Sitter. The former is the conformal or 
strong coupling limit whereas the latter is the decoupling limit. The ghost is absent 
in both cases. Naively, the case f{2H 2 ) = would also appear to be ghost free, 
since the kinetic term for the radion vanishes. Actually, this conclusion is incorrect. 
f(2H 2 ) = corresponds to the case where the radion mixes with the spin 2 mode, 
rendering our analysis invalid. The mixing occurs because the lightest spin 2 mode 
has the same mass as the spin mode {m 2 ight = 2H 2 ). The two modes cease to 
be orthogonal and a more careful analysis is required. This was done for the self- 
accelerating branch of the DGP model, where the ghost was shown to remain even 
when f(2H 2 ) = [10, O, It is natural to expect the same behaviour here. 



6. Discussion 

In this paper we have considered the stability of de Sitter branes in the CGP model: 
an asymmetric generalisation of the DGP model. These vacua include the analogue 
of the normal branch in DGP, as well as the self accelerating branch. Whenever the 
background bulk has infinite volume, we have found, without exception, that linear 
perturbations about these vacua contain ghosts. As for the self accelerating branch of 
DGP, there is always a ghost in either the spin 2 or spin sector. If the spin 2 sector 
is well behaved, there is a spin ghost corresponding to the average radion. If the 
spin sector is well behaved, the helicity-0 part of the lightest spin 2 mode is a ghost. 
A more careful analysis is required in the crossover region, when the two offending 
modes mix with one another. However, our experience from the self-accelerating 
branch of DGP would imply that the ghost remains even in this limit \TU, |I3f| . In 
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the most pathological scenarios, there is yet another ghost corresponding to the 
antisymmetric spin 1 mode. 

It is interesting to note that the only way to avoid ghosts in this model is to 
consider Minkowski branes. This was studied in detail in [E2], where certain inter- 
esting vacua were found to be ghost free. These vacua corresponded to the "stealth" 
models, and had the curious property of giving rise to power law acceleration in the 
presence of matter, before asymptoting to Minkowski space at late times. Indeed, 
the stealth model realises the Cardassian cosmology of Freese and Lewis |29j , as well 



as offering a possible resolution of the coincidence problem. Given these successes of 
these models, it is worth asking whether or not our analysis can shed any light on 
their consistency. 

Of course, the stealth vacua do not include de Sitter branes. In fact, the vacuum 
brane is Minkowski and is known to be ghost free, in contrast to the de Sitter 
branes considered here. What we can say is that the introduction of a small brane 
cosmological constant introduces an instability in the stealth model. It is reasonable 
to extend this conclusion to any type of matter, at least for small H. It follows 
that the stealth model is unstable close to the asymptotically Minkowski limit. The 
question now remains: how dangerous is this instability? 

A ghost will terrorize the vacuum if it couples to ordinary fields. The problem 
is that in a unitary theory, the ghost ought to carry negative energy, and can be 
produced in the vacuum along with ordinary fields without violating energy con- 
servation. In a Lorentz invariant theory, the ghost-non ghost production rate is 
divergent, no matter how weak the coupling! This occurs because one can always 
use Lorentz invariance to perform a boost on the 3-momentum cut-off in loop inte- 
grals. However, a generic Friedmann-Robertson- Walker brane automatically breaks 
Lorentz invariance, so the stealth model does not necessarily suffer from this catas- 
trophic instability (for a related discussion, see f39|). If the ghost only couples weakly 
to other fields, the ghost-non ghost production rate gets suppressed. 

The stealth model contains a decoupling scenario where the would be ghost 
decouples from the spectrum as H — ► 0. This corresponds to the case where we have 
k,L = 0, 0l = — 1 and > 0, 9r = +1, so the cosmological dynamics is governed 
by the following 

p = F(H 2 ) = 6M 2 H 2 - 6M| (^Jh 2 + k 2 R - kj^j + 6M'[H. (6.1) 

For small H, it is easy enough to check that f(2H 2 ) ~ —2M R kn < 0, from which we 
conclude that there is a radion ghost. We know from |22| that the radion decouples 



in the Minkowski limit, so it must be weakly coupled at small H. Given that the 
radion feeds into the brane bending mode, we can see this explicitly by considering 
the coupling of the brane bending mode to matter (see equation ( |4.5| )). The coupling 
strength is given by 1/F'(H 2 ) ~ H/3MI, which does indeed go to zero as H — > 0. 
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We conclude that this particular stealth model will barely be affected by the ghost 
at small H, owing to the weakness of the coupling. At larger values of H, our de 
Sitter brane analysis suggests that the ghost coupling becomes significant, but we 
cannot be sure that these results apply to a general FRW brane. 
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